Optimized Design of Touching Parts of Soil Disinfection Machine Based on Strain Sensing and Discrete Element Simulation

With the increasing level in the intensification of agricultural production in China, continuous cropping obstacles have become a problem that needs to be solved. The use of vertical rotary tillage technology and soil disinfection technology is an effective solution. In this paper, a vertical rotary soil-tilling variable disinfection combine was developed and an on-board control system with STM32 as the control core was designed to realize the real-time acquisition of powder monopoly torque information and the variable application of soil disinfection chemicals. Based on the obtained experimental soil parameters, a discrete element soil particle model was established, and orthogonal experiments were conducted to analyze the single-blade roller tillage process, and the optimal operating parameters were finally selected as 500 mm powder monopoly depth, 320 r/min knife roller speed, and 0.26 m/s forward speed, respectively. The field experiment found that the average tillage depth of the implement was 489 mm, the stability coefficient of tillage depth was 94.50%, the uniformity coefficient of soil disinfection was 85.57%, and the applied amount and the speed ratio coefficient of the given flow were linearly related, respectively. This research provides a technical reference for the deep tillage and soil disinfection of the powder monopoly.


Introduction
Continuous cropping obstacles refer to the abnormal growth and development of crops caused by the continuous cultivation of the same crop or related crops on the same soil, which leads to the formation of a special soil environment and the breeding of a large number of harmful bacteria. As the scale of agricultural production has expanded and the level of intensification has increased in recent years, the problem of continuous crop disorder has become increasingly prominent, which has seriously restricted the sustainable development of the agriculture [1,2]. Since the 1980s, China's vegetable farming industry has developed rapidly. The area of protected vegetable fields has reached 4.8 million hectares, accounting for 30% of the vegetable area, and 3.7% of the national cultivated area, respectively [3]. However, because of continuous cropping for many years, protected vegetable fields have had to endure the breeding and spread of pests and diseases, and the content of the soil organic matter has decreased as a result. Long-term shallow rotary tillage cannot be used to solve the problem of continuous cropping, and, instead, a hard plow bottom layer is formed under the crop root system, and the land becomes less permeable to pests. Furthermore, the air permeability of the land becomes poor, the ability to store water and moisture is reduced, the vegetable growth and development is poor, and the efficiency is significantly reduced [4][5][6][7].
Soil disinfection combined with tillage is an effective solution to continuous cropping obstacles, which can be used to kill most of the virus bacteria in the soil at the same time as cultivation. At present, the vertical spiral blade roll has been widely used in afforestation and reforestation operations due to its reliability and ability to crush soil Based on the above research, it can be inferred that the current soil disinfection machines cannot be used to disinfect and sterilize deep bacteria, and disinfection and sterilization are uneven, while vertical rotary tillage technology can crush the soil deeply, which can promote the penetration of disinfectants and improve the mixing degree of the chemicals and soil, thereby producing a deep uniform disinfection and sterilization. The purpose of this paper was to establish a discrete element soil particle model using EDEM 2018 software, simulate and analyze the process of single blade roll tillage, investigate the influence of different operating parameters on power consumption, and design an onboard control system with STM32 as the control core to achieve the real-time collection of information, such as dust monopoly torque and the variable application of soil disinfection chemicals. The tillage and disinfection performance of the machine was verified through experiments. The objective of this paper was to develop a machine to solve the continuous crop barrier, improve the crop growth environment, and achieve green harvest.

Whole Machine Structure
The structure of the machine designed in this paper is shown in Figure 1, which includes a vertical rotary tillage system and a variable disinfection system. The frame was welded onto the rear of a gearbox. The vertical rotary tillage system is mainly composed of a transmission system and a spiral blade roll. The tillage depth can reach more than 40 cm. The variable disinfection system mainly includes an electric throttle valve, a flow control box, a disinfection spray bar, etc., which together coordinate the operation to achieve a variable disinfection of the soil. disinfection uniformity coefficient reached 64.23%, respectively. Fang et al. designed a solid fumigant disinfector, but their disinfectant could not enter the deeper layers of the soil [30]. Ma et al. used microwave heating of soil to verify its disinfection effect, but no field trials were conducted [31].
Based on the above research, it can be inferred that the current soil disinfection machines cannot be used to disinfect and sterilize deep bacteria, and disinfection and sterilization are uneven, while vertical rotary tillage technology can crush the soil deeply, which can promote the penetration of disinfectants and improve the mixing degree of the chemicals and soil, thereby producing a deep uniform disinfection and sterilization. The purpose of this paper was to establish a discrete element soil particle model using EDEM 2018 software, simulate and analyze the process of single blade roll tillage, investigate the influence of different operating parameters on power consumption, and design an on-board control system with STM32 as the control core to achieve the real-time collection of information, such as dust monopoly torque and the variable application of soil disinfection chemicals. The tillage and disinfection performance of the machine was verified through experiments. The objective of this paper was to develop a machine to solve the continuous crop barrier, improve the crop growth environment, and achieve green harvest.

Whole Machine Structure
The structure of the machine designed in this paper is shown in Figure 1, which includes a vertical rotary tillage system and a variable disinfection system. The frame was welded onto the rear of a gearbox. The vertical rotary tillage system is mainly composed of a transmission system and a spiral blade roll. The tillage depth can reach more than 40 cm. The variable disinfection system mainly includes an electric throttle valve, a flow control box, a disinfection spray bar, etc., which together coordinate the operation to achieve a variable disinfection of the soil.  The machine is connected to the front power machinery through a three-point suspension device, and the angle of the blade into the soil was adjusted by controlling the length of the hydraulic top link. The power output shaft drives the three spur gears located in the same plane behind it through a pair of bevel gears to engage and rotate. The lower side of the spur gear is connected to a spiral blade through a flange plate to realize the vertical rotary tillage of the soil on the silt ridge. The variable disinfection system is mounted onto a rack behind the gearbox. After the spiral blade has been buried, the disinfectant in the disinfection kit can be sprayed onto the soil surface through a filter, reflux pump, and throttle valve in turn. A schematic diagram for soil disinfection is shown in Figure 2.
The machine is connected to the front power machinery through a three-point suspension device, and the angle of the blade into the soil was adjusted by controlling the length of the hydraulic top link. The power output shaft drives the three spur gears located in the same plane behind it through a pair of bevel gears to engage and rotate. The lower side of the spur gear is connected to a spiral blade through a flange plate to realize the vertical rotary tillage of the soil on the silt ridge. The variable disinfection system is mounted onto a rack behind the gearbox. After the spiral blade has been buried, the disinfectant in the disinfection kit can be sprayed onto the soil surface through a filter, reflux pump, and throttle valve in turn. A schematic diagram for soil disinfection is shown in Figure 2.

Design of the Blade Roll Torque Detection System and the Variable Disinfection System
In this study, a variable disinfection system was designed with STM32-M3 as the core, including an upper computer and two lower computers, equipped with a UCOSIII real-time operating system and an EmWin graphical user interface [32]. STC8F2K08S2 was used as the core of the first lower computer to measure the soil cutting torque of the spiral blade roll, while the voltage signal was measured using a full bridge circuit connected by resistive strain gauges, and the torque data were transmitted to the upper computer via a DL-20 wireless communication module. A second lower computer with STM32 as the core was installed in a flow control box, which was equipped with a Hall speed sensor, relay, stepper motor controller, and Hall flow meter, and communication was established with the upper computer through a DL-20 wireless module to complete the real-time regulation and flow control. STM32F103-V3 was used as the core of the upper computer, carrying an EmWin graphical application interface on an LCD screen for the real-time display of tension, flow, speed, and positioning information obtained from a GPS module collected by the lower computer, while the torque information for the machine was visualized with Graph control. In addition, the operator can use a SD card to import the prescription information and continuously extract the disinfection information from the prescription map based on the GPS latitude and longitude and send the control signal to the

Design of the Blade Roll Torque Detection System and the Variable Disinfection System
In this study, a variable disinfection system was designed with STM32-M3 as the core, including an upper computer and two lower computers, equipped with a UCOSIII real-time operating system and an EmWin graphical user interface [32]. STC8F2K08S2 was used as the core of the first lower computer to measure the soil cutting torque of the spiral blade roll, while the voltage signal was measured using a full bridge circuit connected by resistive strain gauges, and the torque data were transmitted to the upper computer via a DL-20 wireless communication module. A second lower computer with STM32 as the core was installed in a flow control box, which was equipped with a Hall speed sensor, relay, stepper motor controller, and Hall flow meter, and communication was established with the upper computer through a DL-20 wireless module to complete the real-time regulation and flow control. STM32F103-V3 was used as the core of the upper computer, carrying an EmWin graphical application interface on an LCD screen for the real-time display of tension, flow, speed, and positioning information obtained from a GPS module collected by the lower computer, while the torque information for the machine was visualized with Graph control. In addition, the operator can use a SD card to import the prescription information and continuously extract the disinfection information from the prescription map based on the GPS latitude and longitude and send the control signal to the second lower computer through a wireless communication mode to realize the variable disinfection operation. This system was also used to set the normal operation flow threshold for an abnormal flow alarm. The overall scheme and system interface of this system are shown in Figures 3 and 4, respectively. second lower computer through a wireless communication mode to realize the variable disinfection operation. This system was also used to set the normal operation flow threshold for an abnormal flow alarm. The overall scheme and system interface of this system are shown in Figures 3 and 4, respectively.

Soil Sample Acquisition
To establish the soil model and ensure the accuracy of the simulation, the soil in the test plot was sampled and then analyzed to obtain the relevant physical properties of the soil, and the characteristic parameters are given in Table 1.
All the soil samples were obtained from the field trial area in Pei County, Xuzhou City, Jiangsu Province, China, and these soil samples were measured using the five-point sampling method. Since this study involved a deep tillage machine, the samples were taken from the surface to 60 cm below it in three layers: upper (0-20 cm), middle (20-40 cm) and lower (40-60 cm), respectively. second lower computer through a wireless communication mode to realize the variable disinfection operation. This system was also used to set the normal operation flow threshold for an abnormal flow alarm. The overall scheme and system interface of this system are shown in Figures 3 and 4, respectively.

Soil Sample Acquisition
To establish the soil model and ensure the accuracy of the simulation, the soil in the test plot was sampled and then analyzed to obtain the relevant physical properties of the soil, and the characteristic parameters are given in Table 1.
All the soil samples were obtained from the field trial area in Pei County, Xuzhou City, Jiangsu Province, China, and these soil samples were measured using the five-point sampling method. Since this study involved a deep tillage machine, the samples were taken from the surface to 60 cm below it in three layers: upper (0-20 cm), middle (20-40 cm) and lower (40-60 cm), respectively.

Soil Sample Acquisition
To establish the soil model and ensure the accuracy of the simulation, the soil in the test plot was sampled and then analyzed to obtain the relevant physical properties of the soil, and the characteristic parameters are given in Table 1. All the soil samples were obtained from the field trial area in Pei County, Xuzhou City, Jiangsu Province, China, and these soil samples were measured using the five-point sampling method. Since this study involved a deep tillage machine, the samples were taken from the surface to 60 cm below it in three layers: upper (0-20 cm), middle (20-40 cm) and lower (40-60 cm), respectively.

Shear Test for the Soil
Samples were obtained from the test field inspection area according to the upper, middle, and lower layers, and four soil samples were taken for each different soil depth. A preloading instrument was used to conduct an 8 h preloading test for the soil samples to ensure that the compactness of the soil samples were the same. The permeable stone was placed into a shear box and then pushed into the soil sample to ensure that its shape was not damaged, followed by the permeable stone and pressure transfer plate and the adjustment of the pressurized frame. A vertical stress of 100, 200, 300, and 400 kPa was applied to each soil layer, the shear rate was set to 0.8 mm/min, and the load was applied uniformly, respectively. The start time at which the force expression starts to change was recorded until either the soil sample was damaged, or the maximum deformation displacement reached 6 mm, and the process of force expression with the displacement change was recorded. The shear strength curves for each soil layer were plotted according to the data collected, as shown in Figure 5.

Depth
Cone Index Samples were obtained from the test field inspection area according to the upper, middle, and lower layers, and four soil samples were taken for each different soil depth. A preloading instrument was used to conduct an 8 h preloading test for the soil samples to ensure that the compactness of the soil samples were the same. The permeable stone was placed into a shear box and then pushed into the soil sample to ensure that its shape was not damaged, followed by the permeable stone and pressure transfer plate and the adjustment of the pressurized frame. A vertical stress of 100, 200, 300, and 400 kPa was applied to each soil layer, the shear rate was set to 0.8 mm/min, and the load was applied uniformly, respectively. The start time at which the force expression starts to change was recorded until either the soil sample was damaged, or the maximum deformation displacement reached 6 mm, and the process of force expression with the displacement change was recorded. The shear strength curves for each soil layer were plotted according to the data collected, as shown in Figure 5. The maximum shear stress for each group of tests was determined using the least squares method to fit the line of the shear strength versus the vertical pressure according to the Coulomb formula, and the shear strength parameters for the different soil layers were finally determined, as shown in Table 2.The internal friction angle is expressed as the slope of the fitted line, and the soil cohesion values displayed in Table 2 is shown as the intersection of the fitted line with the vertical axis. The maximum shear stress for each group of tests was determined using the least squares method to fit the line of the shear strength versus the vertical pressure according to the Coulomb formula, and the shear strength parameters for the different soil layers were finally determined, as shown in Table 2.The internal friction angle is expressed as the slope of the fitted line, and the soil cohesion values displayed in Table 2 is shown as the intersection of the fitted line with the vertical axis.

Soil Strength Test
The prepared soil sample was placed onto the compression plate of the strain-controlled, unconfined compression apparatus, following which the swing wheel was rotated, and the position of the upper compression plate was adjusted to make light contact with the upper surface of the soil sample. The wheel was rotated to keep the strain rate between 1% and 3% per minute, respectively, and the values measured using the dynamometer Sensors 2023, 23, 6369 7 of 21 and the displacement meter were recorded at any time until the values measured by the dynamometer reached a peak value or became constant. Then, a 3% to 5% axial strain was applied to complete this group of tests. The readings obtained from the dynamometer were then used to calculate the pressure for the axial stress using Equation (1): where σ is axial stress, Kpa; C is the dynamometer rate coefficient, N/0.01 mm; R is the dynamometer value, 0.01 mm; and A α is the corrected sample area, cm 2 , respectively. The relationship curve between the axial stress and strain was drawn according to the test calculation results, as shown in Figure 6.
The prepared soil sample was placed onto the compression plate of the strain-controlled, unconfined compression apparatus, following which the swing wheel was rotated, and the position of the upper compression plate was adjusted to make light contact with the upper surface of the soil sample. The wheel was rotated to keep the strain rate between 1% and 3% per minute, respectively, and the values measured using the dynamometer and the displacement meter were recorded at any time until the values measured by the dynamometer reached a peak value or became constant. Then, a 3% to 5% axial strain was applied to complete this group of tests. The readings obtained from the dynamometer were then used to calculate the pressure for the axial stress using Equation (1): where is axial stress, ; is the dynamometer rate coefficient, / 0.01 mm; isthe dynamometer value, 0.01 mm; and is the corrected sample area, cm 2 , respectively.
The relationship curve between the axial stress and strain was drawn according to the test calculation results, as shown in Figure 6.

Soil Particle Modeling
The soil particle shape was equated to a sphere with the radius of 8 mm, and Poisson's ratio was chosen to be 0.35 based on the soil (loam) texture and the empirical equation. The basic parameters for the soil model obtained from the literature and the soil test data are given in Table 3 [29,[33][34][35].

Soil Particle Modeling
The soil particle shape was equated to a sphere with the radius of 8 mm, and Poisson's ratio was chosen to be 0.35 based on the soil (loam) texture and the empirical equation. The basic parameters for the soil model obtained from the literature and the soil test data are given in Table 3 [29,[33][34][35]. The interaction parameters for the soil particles include the static friction coefficient, the rolling friction coefficient, and the recovery coefficient. The determination of this part of the parameters required a large number of experiments to be carried out, and this paper used the soil accumulation angle as the base parameter to determine this part of the interaction parameters. The particle parameters that are similar to those obtained from the soil accumulation angle test were selected in the GMEE material library of EDEM 2018 software, and a virtual simulation model was established to simulate the falling process for the soil in its natural state as an approximate representation of the particle interaction parameters. After several rounds of simulation tests, the interaction parameters for each material were finally determined, as shown in Table 4. According to the above tests, the soil showed a loam texture in the test area, and the soil moisture content in this area was deemed to be high. Therefore, the Hertz-Mindlin with bonding model was selected as the soil particle contact model. In the discrete element software, the model for the bonding bond interaction between these soils includes the following five parameters: normal stiffness, tangential stiffness, critical positive stress, critical shear stress, and bond radius, respectively while the process model for soil particles from the bonded state to failure damage is expressed as follows [36,37]: where F n and F t are the bonding bond normal force, tangential force, N; T n and T t are the bonding bond normal torque, tangential torque, N·m; V n and V t designate the normal velocity, tangential velocity, m/s; ω n and ω t are the normal angular velocity, tangential angular velocity, rad/s; S n and S t represent the normal stiffness, tangential stiffness, N/m 3 ; A is the contact area, m 2 ; J is the moment of inertia, m 4 ; R 1 isthe bonding radius, m; and T is the unit time step, s. Bonding bond breaking condition [36,37]: where σ max -critical normal stress, Pa; and τ max -critical tangential stress, Pa. The bond radius between the soil particles can be calculated from the soil density, according to Equation (4). In addition, the normal stiffness and the critical positive stress were obtained using the unconfined compressive strength test. The critical positive stress was defined as the maximum stress that the soil can withstand in the vertical direction, i.e., σ, when the small principal stress is 0, which was calculated according to Equation (3). Normal stiffness refers to the ability of the vertical direction to resist elastic deformation when subjected to a force and is related to the bond radius of the particles in the discrete element, which was determined using the ratio of the maximum stress to the bond radius in this paper. The critical shear stress refers to the critical positive stress in the soil direct shear test, which is the point at which the soil Coulomb equation is fitted in the unconfined compressive strength test. The tangential stiffness is expressed by the ratio of the maximum shear stress to the bond radius. In summary, the parameters of the Hertz-Mindlin with bonding contact model are shown in Table 5. where ρ s -density of water, g/cm 3 ; ρ w -density of the soil, g/cm 3 ; V s -volume of soil particle water layer, cm 3 ; V w -volume of soil particle, cm 3 ; and R w -radius of a single soil particle, cm, respectively.

Simulation and Postprocessing
According to the operating depth and width of the spiral single blade roll, the soil bin was established, the blade roll model was imported, and the motion parameters of the spiral blade roll were set. The whole simulation process is shown in Figure 7. The simulation shows that the spiral blade roll will form pits in the ground when cutting the soil, and the tilled soil will be scattered in a circular shape across the soil surface with a diameter of approximately 230 mm. As the blade roll moves forward, a large number of soil particles are thrown backwards, filling up the holes again, and as the soil is in a loose state after tillage, and the broken bonding bonds are not recreated, this means that the soil layer is "elevated" as a whole after tillage. Figure 8a shows the cross-sectional view of the soil after plowing, which also confirms that vertical rotary tillage does not break the original soil structure, and helps to maintain the physical and chemical properties of the soil. Figure 8b shows the bounding bond breakage with a constant increment in the number of normal operating orders broken with time, indicating the uniformity of the tillage. The simulation shows that the spiral blade roll will form pits in the ground when cutting the soil, and the tilled soil will be scattered in a circular shape across the soil surface with a diameter of approximately 230 mm. As the blade roll moves forward, a large number of soil particles are thrown backwards, filling up the holes again, and as the soil is in a loose state after tillage, and the broken bonding bonds are not recreated, this means that the soil layer is "elevated" as a whole after tillage. Figure 8a shows the cross-sectional view of the soil after plowing, which also confirms that vertical rotary tillage does not break the original soil structure, and helps to maintain the physical and chemical properties of the soil. Figure 8b shows the bounding bond breakage with a constant increment in the number of normal operating orders broken with time, indicating the uniformity of the tillage.
In this paper, the ratio of the broken bonds to the total bonds in the tillage area was used to represent the crushing rate, and the EDEM postprocessing module 'Set Bin Group' was used to establish the tillage data collection area and divide the collection grid. The bonding bonds formed and broken between all the particles in the collection area were calculated to yield a crushing rate of 92.34%.
The power consumption is an important parameter to measure the performance of the machine [38]. The power consumption of vertical rotary tillage mainly includes three processes: transmission power consumption, forward power consumption, and soil cutting power consumption, and this paper took the spiral single blade roll cutting simulation as the research object to assess the influence of the forward speed, blade roll speed, and tillage depth on the operating power consumption. In the EDEM simulation, analysis of the above two types of power consumption was reflected in the cutting torque and the forward direction resistance suffered by the blade roll. In the postprocessing interface, the total torque applied to the blade roll in the z-direction, and the total resistance applied to the blade roll in the x-direction can be respectively used to generate line graphs according to time steps, as shown in Figures 9 and 10, respectively. cutting the soil, and the tilled soil will be scattered in a circular shape across the soil surface with a diameter of approximately 230 mm. As the blade roll moves forward, a large number of soil particles are thrown backwards, filling up the holes again, and as the soil is in a loose state after tillage, and the broken bonding bonds are not recreated, this means that the soil layer is "elevated" as a whole after tillage. Figure 8a shows the cross-sectional view of the soil after plowing, which also confirms that vertical rotary tillage does not break the original soil structure, and helps to maintain the physical and chemical properties of the soil. Figure 8b shows the bounding bond breakage with a constant increment in the number of normal operating orders broken with time, indicating the uniformity of the tillage. In this paper, the ratio of the broken bonds to the total bonds in the tillage area was used to represent the crushing rate, and the EDEM postprocessing module 'Set Bin Group' was used to establish the tillage data collection area and divide the collection grid. The bonding bonds formed and broken between all the particles in the collection area were calculated to yield a crushing rate of 92.34%.
The power consumption is an important parameter to measure the performance of the machine [38]. The power consumption of vertical rotary tillage mainly includes three  processes: transmission power consumption, forward power consumption, and soil c ting power consumption, and this paper took the spiral single blade roll cutting simu tion as the research object to assess the influence of the forward speed, blade roll spe and tillage depth on the operating power consumption. In the EDEM simulation, analy of the above two types of power consumption was reflected in the cutting torque and t forward direction resistance suffered by the blade roll. In the postprocessing interface, t total torque applied to the blade roll in the z-direction, and the total resistance applied the blade roll in the x-direction can be respectively used to generate line graphs accordi to time steps, as shown in Figures 9 and 10, respectively. The above figure shows that in the time period of 0~2 s particle bed generation, t torque and forward resistance are 0; in the time period of 2~4 s, the spiral blade rolls in the soil, the torque increases rapidly, and the tension swings back and forth between ±5 N; and in the time period of 4~10 s, the normal operation time for the blade roll, the torq gradually stabilizes over a certain interval, while the tension swings rapidly to the m significant value point due to the beginning of the forward movement, and remains sta over this time interval.  The above figure shows that in the time period of 0~2 s particle bed generation, the torque and forward resistance are 0; in the time period of 2~4 s, the spiral blade rolls into the soil, the torque increases rapidly, and the tension swings back and forth between ±500 N; and in the time period of 4~10 s, the normal operation time for the blade roll, the torque gradually stabilizes over a certain interval, while the tension swings rapidly to the most significant value point due to the beginning of the forward movement, and remains stable over this time interval.

Validation of the Simulation System
To ensure the accuracy of the simulation results, the torque was assessed using the resistance-strain type measurement method, while the torque value was calculated using the changing resistance value, and the calculated results were compared with the simulation results to verify whether this simulation system is reliable. The actual operating condition was integrated, and a position 30 mm below the reinforced rib plate was selected as the patching site. The small axis calibration method was used to calibrate the pasted strain gauge torque sensor. By hanging weights onto the end of the arm bar, a torque was generated on the blade roll, and the strain voltage output from the bridge was measured using an oscilloscope. The calibration results are shown in Figure 11.

Validation of the Simulation System
To ensure the accuracy of the simulation results, the torque was assessed using the resistance-strain type measurement method, while the torque value was calculated using the changing resistance value, and the calculated results were compared with the simulation results to verify whether this simulation system is reliable. The actual operating condition was integrated, and a position 30 mm below the reinforced rib plate was selected as the patching site. The small axis calibration method was used to calibrate the pasted strain gauge torque sensor. By hanging weights onto the end of the arm bar, a torque was generated on the blade roll, and the strain voltage output from the bridge was measured using an oscilloscope. The calibration results are shown in Figure 11.

Validation of the Simulation System
To ensure the accuracy of the simulation results, the torque was assessed using the resistance-strain type measurement method, while the torque value was calculated using the changing resistance value, and the calculated results were compared with the simulation results to verify whether this simulation system is reliable. The actual operating condition was integrated, and a position 30 mm below the reinforced rib plate was selected as the patching site. The small axis calibration method was used to calibrate the pasted strain gauge torque sensor. By hanging weights onto the end of the arm bar, a torque was generated on the blade roll, and the strain voltage output from the bridge was measured using an oscilloscope. The calibration results are shown in Figure 11.  The torsion measuring device was fixed with tape, and the data were transmitted to a PC via a serial port. A single blade roll was installed onto the machine, and the torsion measuring device was then added. The simulation operation process was simulated using the motion parameters at the time of the simulation to test the single blade roll in the field, and the operation process is shown in Figure 12. The torsion measuring device was fixed with tape, and the data were transmitted to a PC via a serial port. A single blade roll was installed onto the machine, and the torsion measuring device was then added. The simulation operation process was simulated using the motion parameters at the time of the simulation to test the single blade roll in the field, and the operation process is shown in Figure 12. To cope with the sudden influx of data with positive and negative oscillations and increase the data processing time, the reference voltage was set at 12 mV, the expected bit value was 768, and the actual bit value was at 771, respectively. The measured voltage signal bit value graph is shown in Figure 13. The actual working state bit average value was 2381, and the difference was 1610, corresponding to a voltage signal of 25.15625 mV and a torque of 157.41, respectively. The simulated torque value was 148.02 N × m, and the error value was 5.97%, respectively. The reason for this analysis was that the field conditions are more complex, and the weeds and debris among the soil and rocks can produce errors in the experimental results. In summary, these experimental results are basically consistent with the simulation results, thereby indicating that the simulation model established in this study is feasible. To cope with the sudden influx of data with positive and negative oscillations and increase the data processing time, the reference voltage was set at 12 mV, the expected bit value was 768, and the actual bit value was at 771, respectively. The measured voltage signal bit value graph is shown in Figure 13. The actual working state bit average value was 2381, and the difference was 1610, corresponding to a voltage signal of 25.15625 mV and a torque of 157.41, respectively. The simulated torque value was 148.02 N × m, and the error value was 5.97%, respectively. The reason for this analysis was that the field conditions are more complex, and the weeds and debris among the soil and rocks can produce errors in the experimental results. In summary, these experimental results are basically consistent with the simulation results, thereby indicating that the simulation model established in this study is feasible.

Virtual Simulation Experiment
The orthogonal test is an effective method for constructing multi-factor and multilevel tests [39]. From the perspective of practical simulation research, the forward speed,

Virtual Simulation Experiment
The orthogonal test is an effective method for constructing multi-factor and multi-level tests [39]. From the perspective of practical simulation research, the forward speed, the blade roll speed, and the tillage depth were selected as the test factors, and the power consumption for soil cutting and the forward power consumption were taken as the test indices. According to the pre-simulation test and the vertical rotary tillage requirements, with the actual operation standards, the test factor level table was set, as shown in Table 6. In this paper, a single blade roll was removed to simulate the operation process and the magnitude of the power consumption of the computer tool. This virtual simulation test was used to mainly study the forward power consumption and the soil cutting power consumption of the spiral blade roll in the operating state of vertical rotary tillage, and the soil cutting power consumption was calculated according to Equation (5): where P z -cutting power consumption, kW; T z -cutting torque, N·m; and n 1 -blade roll speed, r/min, respectively. Forward power consumption is the rate of energy consumed by the machine to keep moving in the forward direction, and was calculated using Equation (6): where P x -forward power consumption, W; F x -rotary tillage pull, N; and V x -forward speed, m/s, respectively. The three-factor, three-level orthogonal virtual test was designed according to Table 6 [40,41], where A represents the forward speed; B represents the tillage depth; and C represents the blade roll speed, respectively. The test results are shown in Table 7. It can be observed from the analysis table that within the range of the set test group, the total power consumption for single blade roll vertical rotary tillage basically ranges between 2500-6500 W, while the vertical rotary soil-tilling variable disinfection combine designed in this paper is composed of three blade rolls, and its rotary tillage power consumption varies within 20 kW. According to the above analysis, the proportion of the forward power consumption to the total power consumption of the ridge was 4.5~12.5%, which is basically consistent with the proportion of the power consumption for the traditional rotary cultivators.

Analysis of the Range
To analyze the role of each factor on the power consumption, the power consumption for soil cutting along with the forward power consumption was determined with range analysis through the R value obtained from the size of the analysis of the influence in the order of priority, and the range analysis table is shown in Table 8. According to the range analysis for the two types of power consumption, the factors that affect the power consumption performance of the machine follow the order of tillage depth, forward speed, and blade roll speed, respectively. Most of the power consumption for the tillage implements arise from the soil cutting, so soil cutting power consumption is the focus of the research in the study of tillage power consumption. To further investigate the influence of each operating parameter on the power consumption for soil cutting, a three-factor ANOVA was conducted for the power consumption of soil cutting using SPSS 22.0 software with the three factors of forward speed, tillage depth, and blade roll speed as variables, and Table 9 shows the ANOVA table obtained for the three-factor test. Table 9. Three-factor ANOVA (analysis of variance). The ANOVA results are consistent with the results of the ANORA. The p values for the forward speed and tillage depth are 0.029 and 0.012, respectively, which indicates that the forward speed and the tillage depth have significant effects on the power consumption. The power consumption for soil cutting, the operation quality, the blade roll operation intensity, and the agronomic requirements for vertical rotary tillage were integrated; the selected tillage depth was 50 cm, the forward speed was 0.26 m/s, and the blade roll speed was not determined.

Single-Factor Test for the Blade Roll Speed
To investigate the effects of the blade roll speed on the power consumption for soil cutting, single-factor tests were designed for analysis. Five sets of single-factor tests were designed with a fixed forward speed of 0.26 m/s and a tillage depth of 500 mm in the speed range of 260-340 r/min, respectively. The experimental data were calculated to derive the average soil cutting power consumption and the torque variance. The changes in the soil cutting power consumption and the operational stability with a change in the blade roll speed were analyzed, and the results are shown in Figure 14.

Tillage Performance Field Test
For a vertical rotary tillage machine, the tillage depth and the tillage stability are the most important parameters. The length of the stable working area was set to 12 m, the working depth was measured by using a steel ruler, and 29 sample points were taken from the stable working area in a serpentine way to measure the tillage depth. The actual point-taking diagram is shown in Figure 15.  The analysis shows that the power consumption for soil cutting rises slowly with increasing rotational speed, but the trend was not determined to be significant, which verifies the accuracy of the above orthogonal test. The torque variance reflects the stability of the machine operation and the impact load problem. The smaller the variance, the more stable the operation of the implement, and the smaller the impact load on it. Therefore, the final combination of optimal tillage parameters was determined to be 500 mm tillage depth, 0.26 m/s forward speed, and 320 r/min blade roll speed, respectively.

Tillage Performance Field Test
For a vertical rotary tillage machine, the tillage depth and the tillage stability are the most important parameters. The length of the stable working area was set to 12 m, the working depth was measured by using a steel ruler, and 29 sample points were taken from the stable working area in a serpentine way to measure the tillage depth. The actual point-taking diagram is shown in Figure 15.
For a vertical rotary tillage machine, the tillage depth and the tillage stability are the most important parameters. The length of the stable working area was set to 12 m, the working depth was measured by using a steel ruler, and 29 sample points were taken from the stable working area in a serpentine way to measure the tillage depth. The actual point-taking diagram is shown in Figure 15. The average tillage depth, the standard deviation of the tillage depth, the coefficient of variation of the tillage depth, and the stability coefficient of the tillage depth for the machine were all calculated using the following equations: where ̅ is the average tillage depth, mm; is the tillage depth of the th sample, mm; and represents the sample size, 29, respectively. Equation (8) was used to determine the standard deviation in the tillage depth. The average tillage depth, the standard deviation of the tillage depth, the coefficient of variation of the tillage depth, and the stability coefficient of the tillage depth for the machine were all calculated using the following equations: where S is the average tillage depth, mm; S i is the tillage depth of the ith sample, mm; and N represents the sample size, 29, respectively. Equation (8) was used to determine the standard deviation in the tillage depth.
Equation (9) was used to determine the coefficient of variation for the tillage depth.
Equation (10) was used to determine the stability coefficient for the tillage depth.

Variable Disinfection Field Test
Eight hundred grams of high potassium water-soluble fertilizer was dissolved in 25 L of water and mixed well as a disinfection agent for spraying, and the variable disinfection effect was assessed by detecting the content and distribution of potassium elements in the soil before and after machine operation. Before the operation, the soil area to be tested was sampled according to the serpentine sampling points to detect the amount of potassium originally contained in the soil to be sprayed, and the samples were taken according to the surface soil and subsoil, with the surface soil obtained directly from the ground, and the subsoil sample taken at a distance of 15 cm below the ground surface. A flow rate ratio of four was set, and the test area was then operated on. The soil samples were resampled after 12 h of resting to form a control with the previous sampling points, and the potassium content of the soil samples was measured using a soil nutrient quick tester after the sampling was completed. The potassium content for each point was obtained using a soil nutrient tester, and the results are shown in Figure 16. Analysis of the actual potassium application to the soil shows that the potassium water-soluble fertilizer can easily penetrate into the lower soil layer. The uniformity coefficient and the application amount distribution stability coefficient were calculated for the application amount of the disinfecting agent to the surface layer of the soil, as according to Equations (11)-(13).
= (100 − ) × 100% Substituting the actual application rates for the above topsoil into the above equation results in a standard deviation of 15.99 for the distribution of applied potassium in the topsoil, a coefficient of variation of distribution V = 14.43, and a coefficient for the stable nature of the distribution of the application rate of U = 85.57%, respectively.
To further evaluate the accuracy of the spraying flow rate, the flow rate was adjusted to a flow rate ratio of 3 and 5, respectively. The above process was repeated, and only the Analysis of the actual potassium application to the soil shows that the potassium water-soluble fertilizer can easily penetrate into the lower soil layer. The uniformity coefficient and the application amount distribution stability coefficient were calculated for the application amount of the disinfecting agent to the surface layer of the soil, as according to Equations (11)- (13).
Substituting the actual application rates for the above topsoil into the above equation results in a standard deviation of 15.99 for the distribution of applied potassium in the topsoil, a coefficient of variation of distribution V = 14.43, and a coefficient for the stable nature of the distribution of the application rate of U = 85.57%, respectively.
To further evaluate the accuracy of the spraying flow rate, the flow rate was adjusted to a flow rate ratio of 3 and 5, respectively. The above process was repeated, and only the amount of potassium applied to the topsoil was assessed; the results are shown in Table 11. Table 11. Potassium application at different rate ratios.

Flow Rate Ratio of 3 Flow Rate Ratio of 5 Sample Points Potassium Application/ppm Sample Points Potassium Application/ppm
The average value for the flow rate ratio was extracted, and the three sets of test dat with coefficients of three, four, and five were removed. The data point with a flow rat ratio of four was used as the basis to draw the theoretical application line, and the posi tions of the points with flow rate ratios of three and five were observed, as shown in Figur 17.

Discussion
In this research, we aimed to develop a soil disinfection machine and assess its ra tionality by the machine's tillage performance and the disinfection effect observed in th

Discussion
In this research, we aimed to develop a soil disinfection machine and assess its rationality by the machine's tillage performance and the disinfection effect observed in the field. The results of the simulation and field trial showed that the machine could achieve a tillage depth stability factor of 94.50% and a disinfection uniformity factor of 85.57%, respectively. Compared with the vertical spiral blade roll used in this research, Zhang and Liu et al. used multiple pairs of square knives as rototiller rolls to tillage the soil [42], and their tillage depth stability reached 95.24% and 94.84%, respectively, which were close to the results of this research. However, since this type of blade roll throws soil in all directions, which can affect the mixing of the disinfection chemicals with the soil, and with their disinfection uniformity coefficient being 64.23%, it can be shown that the vertical spiral blade roll is more suitable for vertical tillage as well as soil disinfection.

Conclusions
In this paper, a vertical rotary tilling-soil variable disinfection combine was designed, simulated, and tested. The study shows that this machine is capable of vertical rotary tillage operation while achieving the variable chemical disinfection agent spraying of the soil. After determining the soil parameters, an accurate discrete element soil particle model was established. A spiral blade roll-soil interaction model was established based on the discrete element method, and the accuracy of the prototype was verified through applying field tests. The results show that the relative error between the actual soil cutting torque and the simulation results under the same motion parameters was 5.97%. The simulation system was also used to conduct orthogonal tests, and the optimal operating parameters were finally selected from the viewpoint of reducing the operating power consumption and improving the operating effect: tillage depth of 500 mm, blade roll speed of 320 r/min, and a forward speed of 0.26 m/s, respectively. A blade roll torque detection system and a variable disinfection system based on STM32 was designed to achieve a variable disinfection operation for the soil through the synergistic cooperation of the GPS and a prescription map with an electronically controlled throttle valve as the main actuating device. The test results show that the uniformity coefficient for disinfection in the upper soil layer reached 85.57%, the disinfection agent can penetrate to a distance of 15 cm below the surface, and the application amount and given flow rate ratio coefficients are basically linear.
In this study, the design of the spiral blade roll was only based on the empirical formula and the existing research base. The parameters, such as the spiral lift angle of the blade roll have an important influence on the power consumption of the machine, as well as the operation quality. In future research, we will use discrete element software to perform soil cutting simulations for different parameters of the spiral knife rolls, including variable pitch, segmented spiral, and double spiral blade rolls to achieve further optimization of the implements. In addition, the variable disinfection method used in this study does not form a closed-loop control, although a subsequent information collection module such as a flow meter was added. In future research, we will establish a PID closed-loop control model to correct the actual output flow rate by imposing a corrective control relationship using the flowmeter signal as feedback information.